Method and apparatus for performing motion compensated temporal filtering in video encoding

ABSTRACT

A method ( 50 ) and apparatus ( 40 ) for reducing noise in a video signal calculates ( 53 ) a discrimination value of a filter based on a magnitude difference between a current frame and a prior frame of a sequence of frames and modifies ( 54 ) the current frame using the modified discrimination value. The output of the filter is calculated as follows:
 
 g ( i, j, k )= f ( i, j, k )±β
 
wherein:
         f (i, j, k) is the current frame input to the filter;   g(i, j, k) is the current frame output from the filter; and   β is the discrimination value of the filter, which is calculated as follows:       

     
       
         
           
             β 
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           κδ 
                           ′ 
                         
                         + 
                         χ 
                       
                       , 
                     
                   
                   
                     
                       
                         
                           δ 
                           ′ 
                         
                         = 
                         
                           max 
                           ⁢ 
                           
                             { 
                             
                               δ 
                               , 
                               
                                 δ 
                                 0 
                               
                             
                             } 
                           
                         
                       
                       , 
                     
                   
                   
                     
                       δ 
                       &lt; 
                       
                         δ 
                         1 
                       
                     
                   
                 
                 
                   
                     
                       0 
                       , 
                     
                   
                   
                     
                         
                     
                   
                   
                     
                       δ 
                       ≥ 
                       
                         δ 
                         1 
                       
                     
                   
                 
               
             
           
         
       
         
         
           
             
               
                 in which: 
                 κ and χ are constants less than one; 
                 δ is a magnitude difference between the current frame and the prior frame; and 
                 δ 0  and δ 1  are constants that satisfy δ 0 ≦δ 1 .

FIELD OF THE INVENTION

The present invention relates generally to methods and apparatuses forprocessing digital information, and more particularly to a method andapparatus for processing digital information prior to transmission overa network.

BACKGROUND

Motion compensated temporal filters (MCTF) are used in pixel processingof video images. MCTF is an adaptive temporal filter along estimatedmotion trajectories. The temporal correlation is exploited throughtemporal filtering. MCTF is most effective in processing sequences ofimages in which signal is highly temporal correlated along the motiontrajectories but the noise, by nature, is totally uncorrelated. In otherwords, the image does not change from frame to frame but the noise does.The noise probability distribution is also exploited. MCTF is primarilyused to reduce the zero mean stationary white Gaussian noise. Such noisehas a bell shaped distribution, and it is of primarily small magnitude.Large magnitude noise rarely happens. Based on these two observations,the magnitude differences between the current frame and the previousfiltered frame are first calculated along the estimated motiontrajectories. Small differences are considered primarily of noise. Largedifferences are considered primarily of signal content.

However, existing MCTFs do not remove all of the noise inherent in avideo signal, thereby wasting valuable bandwidth resources.

The present invention is therefore directed to the problem of developinga method and apparatus for performing filtering of noise in a videoencoding process that removes as much of the noise as possible withoutdisrupting the video image.

SUMMARY OF THE INVENTION

The present invention solves these and other problems by providing amethod and apparatus for filtering noise in a video encoding process, inwhich an adaptation scheme based on differences between frames isemployed, thereby allowing the signal content to pass through whilestill effectively suppressing the random noise.

According to one aspect of the present invention, a method for filteringnoise in a video encoding process determines a discrimination value of afilter based on a magnitude difference between a current frame and aprior frame of a sequence of frames and modifies the current frame usingthe determined discrimination value.

According to another aspect of the present invention, the current framemay be modified as follows:g(i, j, k)=f(i, j, k)±βwherein f(i, j, k) is the current frame input to the filter; g(i, j, k)is the current frame output from the filter; and β is the discriminationvalue of the filter.

According to still another aspect of the present invention, thediscrimination value may be determined as follows:

$\beta = \{ \begin{matrix}{{\delta - {( {1\text{/}\Delta} )\delta^{2}}},} & {\delta < \Delta} \\{0,} & {\delta \geq \Delta}\end{matrix} $in which δ is the magnitude difference between the current frame and theprior filtered frame along the motion trajectories, and Δ is apredetermined filter strength.

According to yet another aspect of the present invention, thediscrimination value may be determined as follows:

$\beta = \{ \begin{matrix}{{{\kappa\delta}^{\prime} + \chi},} & {{\delta^{\prime} = {\max\{ {\delta,\delta_{0}} \}}},} & {\delta < \delta_{1}} \\{0,} & \; & {\delta \geq \delta_{1}}\end{matrix} $in which κ and χ are constants less than one; δ is a magnitudedifference between the current frame and the prior filtered frame alongthe motion trajectories; δ₀ and δ₁ are constants that satisfy δ₀≦δ₁; andκ, χ, δ₀, and δ₁ can be modified to adjust the filter strength.

According to still another aspect of the present invention, an apparatusfor encoding a video signal includes two encoders and a noise reductionpreprocessor. The first encoder receives the video signal. The secondencoder receives a plurality of parameters from the first encoder to usein encoding the video signal. The noise reduction preprocessor receivespixel data from the first encoder, filters the pixel data by determininga discrimination value of a filter based on a magnitude differencebetween a current frame and a prior frame of a sequence of frames, andby modifying the current frame using the modified discrimination value,and then provides preprocessed video to the second encoder.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts a plot of a linear adaptive filter used in a MotionCompensation Temporal Filter (MCTF) and the optimal filter according toone aspect of the present invention.

FIG. 2 depicts a mapping of the linear adaptive filter in FIG. 1 to adifferent domain and a plot of an exemplary embodiment of a MCTFaccording to another aspect of the present invention.

FIG. 3 depicts the exemplary embodiment of FIG. 2 in the original domainalong with the linear adaptive filter and the optimal filter.

FIG. 4 depicts a block diagram of an exemplary embodiment of a pixelprocessing system, which employs the MCTF of the present inventionaccording to still another aspect of the present invention.

FIG. 5 depicts a flow chart of an exemplary embodiment of a method forfiltering images in a video signal according to yet another aspect ofthe present invention.

FIG. 6 depicts a block diagram of an exemplary embodiment of a MCTF forfiltering images in a video signal according to yet another aspect ofthe present invention.

FIG. 7 depicts a plot of a processing order used in a bi-directionalMCTF according to yet another aspect of the present invention.

DETAILED DESCRIPTION

It is worthy to note that any reference herein to “one embodiment” or“an embodiment” means that a particular feature, structure, orcharacteristic described in connection with the embodiment is includedin at least one embodiment of the invention. The appearances of thephrase “in one embodiment” in various places in the specification arenot necessarily all referring to the same embodiment.

MCTF is used in a dual-pass video encoder. MCTF is employed between thefirst pass and second pass encoders. The first pass encoder providesmotion vectors, film mode status, scene change decisions, andcomplexities along with the video data to MCTF. MCTF filters the videodata using the motion vectors, film mode status, scene change decisions,and complexities from first pass encoder. MCTF then provides thefiltered video data to second pass encoder, which also receivesinformation from first pass encoder.

FIG. 4 shows the block diagram of the overall encoding process forperforming MPEG2 encoding. The architecture implements certain videocompression features to improve the video quality over current videoencoder products, including motion compensated temporal filtering thattracks motion of objects to efficiently remove video noise. There arefive major components in the system—a first pass encoder 41, a noisereduction preprocessor 42, a second pass encoder 43, a mastercompression controller 44 and a packet processor 45. There is also aVideo FIFO and Packet Creator 46.

The first pass encoder 41 is a 110 MHz C-Cube DvXpertII (E4) MPEG2 videocompression chip. It performs the following functions:

-   -   Captures the input CCIR656 (D1) format video;    -   Captures VBI data;    -   Horizontal filters and decimates the pixel data in ½, ⅔, or ¾        resolution mode;    -   Detelecines (inverse 3:2 pulldown);    -   Detects special events (scene changes);    -   Trial encodes the video frames to generate complexity measure        for the frames; and    -   Determines P or B picture coding type of a frame.

The noise reduction preprocessor 42 is a 300 MHz Equator MAPCA mediaprocessor chip. It performs the following functions:

-   -   Lookahead frame delay;    -   Temporal filtering; and.    -   Spatial filtering.

The motion compensated temporal filter is an IIR filter. The motionvectors generated by first pass encoding are used to trace motion ofobjects.

FIG. 6 shows the block diagram of the MCTF. The forward temporal filter61 processes each frame (P and B) in the video sequence using theforward motion vectors from the first pass encoder 41. After a certaindelay, the backward temporal filter 63 processes the B frames in thevideo sequence using the backward motion vectors from the first passencoder 41. If there are two B frames between two consecutive P frames,minimum two-frame delay is required between the forward 61 and backward63 temporal filters. The backward temporal filter 63 does not processthe P frames in the video sequence.

For P frames, unidirectional temporal filtering is employed, the forwardtemporal filter generates the following output g_(P)(, j, k) for everypixel at indexes (i, j) on (x, y) coordinate and index k on time t axis.g _(P)(i, j, k)=f _(p)(i, j, k)±β

-   -   where        -   f_(p)(i, j, k) is the input pixel of the current P frame at            input time k.        -   β is the discrimination value of the filter and it is a            function of the absolute difference between f_(P)(i, j, k)            and g_(P)(i+d_((x, f)), j+d_((y, f)), k−M−1).        -   (d_((x, f)), d_((y, f)) is the forward motion vector at            (i, j) with half-pel precision.            g_(P)(i+d_((x, f)), j+d_((y, f)), k−M−1) is the prior            filtered P frame along the motion trajectories at input time            k−M−1.        -   M is the number of B frames between the two consecutive P            frames.

For B frames, bi-directional temporal filtering is employed, the forwardtemporal filter generates the following output g′_(B)(i, j, k) for everypixel at indexes (i, j) on (x, y) coordinate and index k on time t axis.g′ _(B)(i, j, k)=f _(B)(i, j, k)±β₁

-   -   where        -   f_(B)(i, j, k) is the input pixel of the current B frame at            input time k.

For the first B frame between the two consecutive P frames, β₁ is thediscrimination value of the filter and it is a function of the absolutedifference between f_(B)(i, j, k) and g_(P)(i+d_((x, f)),j+d_((y, f) k−)2).

For the second B frame between the two consecutive P frames, β₁ is thediscrimination value of the filter and it is a function of the absolutedifference between f_(B)(i, j, k) and g_(P)(i+d_((x, f)), j+d_((y, j)),k−2).

(d_((x, f)), d_((y, f))) is the forward motion vector at (i, j) withhalf-pel precision. g_(P)(i+d_((x, f)), j+d_((y, f)), k−1) org_(P)(i+d_((x, f)), j+d_((y, f)), k−2) is the prior forward temporalfiltered P frame along the motion trajectories at input time k−1 or k−2respectively.

For B frames, the backward temporal filter generates the final outputg_(B)(i, j, k) for every pixel at indexes (i, j) on (x, y) coordinateand index k on time t axis.g _(B)(i, j, k)=g′ _(B)(i, j, k)±β₂

-   -   where        -   g′_(B)(i, j, k) is the input pixel to the backward temporal            filter of the current B frame at input time k. g′_(B)(i,            j, k) has been previously filtered by the forward temporal            filter.

For the first B frame, β₂ is the discrimination value of the filter andit is a function of the absolute difference between g′_(B)(i, j, k) andg_(P)(i+d_((x, b)), j+d_((y, b)), k+2).

For the second B frame, β₂ is the discrimination value of the filter andit is a function of the absolute difference between g′_(B)(i, j, k) andg_(P)(i+d_((x, b)), j+d_((y, b)), k+1).

(d_((x, b)), d_((y, b))) is the backward field motion vector at (i, j)with half-pel precision.

g_(P)(i+d_((x, b)), j+d_((y, b)), k+1) or g_(P)(i+d_((x, b)),j+d_((y, b)), k+2) is the prior forward temporal filtered P frame alongthe motion trajectories at input time k+1 or k+2 respectively.

It is important to insert a minimum two-frame delay between the forwardand backward temporal filters. The two-frame delay allows the forwardtemporal filter 61 process the P frame at input time k+1 or k+2 prior tothe backward temporal filter 63 process the B frame at input time k.

FIG. 7 shows the order of processing for the bi-directional filtering.The forward temporal filter processes each frame (both P and B) in thedisplay order as shown by the arrows 1, 2, 3, and 4. The backwardtemporal filter processes the B frames in the display order as shown bythe arrows 5 and 6. It does not process the P frames.

The second pass encoder 43 is a 110 MHz C-Cube DvXpertII (E4) MPEG2video compression chip. It performs the following functions:

-   -   Refines the motion vectors generated by the first pass encoding;    -   Calculates need parameters for statmux;    -   Performs MPEG2 Main Profile at Main Level video encoding;    -   Rate control; and    -   Determine if a frame should be coded as I-frame.

The Master Compression Controller (MCC) 44 is a 50 MHz PPC860T PowerPCcontroller. It performs the following functions.

-   -   Communicates with the Controller/Element Manager;    -   Functions as the host processor for the first pass encoder,        second pass encoder and noise reduction preprocessor to perform        configuration and status inquiry tasks;    -   Decodes the VBI data captured on the first pass encoder, delays        the VBI decoded VBI data and passes it to the second pass        encoder;    -   Communicates with the QLP to relay need parameter and bit rate        messages to the second pass encode; and    -   Maintains the queue of bit rate values received from the QLP for        buffer modeling.

The Packet Processor is an 80 MHz DSP56309. It performs the followingfunctions.

-   -   Generates PTS and DTS;    -   Pulls bits out from the video FIFO to create MPEG2 packets at a        fixed rate, or at a bit rate specified by the transmission bit        rate parameter received from the QLP;    -   Models the decoder buffer and computes the minimum and maximum        transmission bit rate constraints to send to the QLP in statmux;        and    -   Performs other MPEG2 system multiplex functions.

According to one aspect of the present invention, half-pel motioncompensation is employed in the MCTF to improve the accuracy of theestimated motion trajectories. This increases the signal correlation toallow more noise reduction without trailing or blocking artifacts.

According to another aspect of the present invention, the adaptationscheme of the MCTF is optimized. The optimal adaptation scheme is firstdetermined by experiment. By adjusting the filter continuously atvarious levels using a series of images with variety of contents, theoptimal adaptation area is identified. The “forbidden area” shown inFIG. 1 represents the area in which filtering would degrade the image.The area in white represents the area in which potential noise that maybe filtered without adversely affecting the image. The optimaladaptation allows more noise reduction without adversely affecting thesignal content. In one version of MCTF, a linear adaptation scheme isemployed. Another MCTF is developed, which employs a non-linearadaptation scheme to closely approximate the optimal scheme, and whichalso reduces the computational complexity if parameters are selectedappropriately.

The Optimal Adaptation Scheme

The MCTF is described by the following equation:g(i, j, k)=αf (i, j, k)+(1−a)g(i−d _(x) , j−d _(y) , k−1)  (1)where g(i−d_(x), j−d_(y), k−1) is the previous filtered frame; f(i, j,k) and g(i, j, k) are the current input and output frames, respectively;d_(x) and d_(y) are the horizontal and vertical motion displacements,respectively; the filter coefficient α adapts to the magnitudedifference δ between f(i, j, k) and g(i−d_(x), j−d_(y), k−1). Theoptimal adaptation scheme is established experimentally to allow maximumnoise suppression without adversely affecting the video content. Theoptimal adaptation scheme is shown in FIG. 1 in the α-δ domain alongwith the maximum strength linear adaptation scheme. From FIG. 1, we cansee the optimal adaptation scheme is a non-linear function which dependson δ^(γ), where γ<1. Obviously, the implementation of the optimaladaptation scheme is computationally intensive in the α-δ, domainbecause of the need to calculate δ^(γ). Also, any adaptation schemecannot travel through the shaded area. Otherwise, it will adverselyaffect the video content.

While computationally intensive, the non-linear adaptation scheme willresult in an improvement in noise reduction. Nevertheless, the followingexemplary embodiment seeks to obtain this improvement whileconcomitantly decreasing the computation complexity.

Exemplary Embodiment

An exemplary embodiment of a Motion Compensated Temporal Filter (MCTF)is developed in this section. Let δ represent the magnitude differencebetween f(i, j, k) and g(i−d_(x), j−d_(y), k−1).δ=|f(i, j, k)−g(i−d _(x) , j−d _(y) , k−1)| or g(i−d _(x) , j−d _(y) ,k−1)=f(i, j, k)±δ  (2)

Substituting equation (2) into equation (1), we obtain:g(i, j, k)=f(i, j, k)±(1−α)δ=f(i, j, k)  (3)where β represents the discrimination between the filtered and originalframes.β=(1−α)δ  (4)

Equation (3) represents a new implementation of a MCTF, where thediscrimination β is first calculated. Adding or subtracting β from theinput, we can then obtain the filter output. If f(i, j, k) is greaterthen g(i−d_(x), j−d_(y), k−1), subtraction is performed. Otherwiseaddition is performed. Comparing equation (3) to equation (1), we cansee equation (1) consists of two multiplications while equation (3)consists of only one. Because α is a fixed-point number and the pixelvalues are 8-bit precision, the multiplications have to be performedwith 16-bit precision. Although equation (3) has only onemultiplication, the additional computation involved in the decision ofwhether adding or subtracting β from the input cancels out the savings.

Equation (4) translates an adaptation scheme in the original α-δ domaininto the new β-δ domain. As will be shown, a non-linear adaptation inthe α-δ domain is translated into a linear adaptation in the β-δ domain.This allows a non-linear adaptation in the original domain to beimplemented as a linear adaptation in the new domain. This is thefundamental advantage of the exemplary embodiment.

To get more insight into the exemplary embodiment, let us first considerthe linear adaptation scheme depicted in FIG. 1. In the linearadaptation scheme, α is a linear function of δ.

$\begin{matrix}{\alpha = \{ \begin{matrix}{{( {1\text{/}\Delta} )\delta},} & {\delta < \Delta} \\{0,} & {\delta \geq \Delta}\end{matrix} } & (5)\end{matrix}$

Clearly, the bigger is Δ, the stronger is the filter. Substituteequation (5) into equation (4), we obtain:

$\begin{matrix}{\beta = \{ \begin{matrix}{{\delta - {( {1\text{/}\Delta} )\delta^{2}}},} & {\delta < \Delta} \\{0,} & {\delta \geq \Delta}\end{matrix} } & (6)\end{matrix}$

Equation (6) describes the adaptation scheme in the new β-δ domain.Clearly, it is a non-linear function. In another word, a linearadaptation scheme in the α-δ domain (original implementation) istranslated into a non-linear adaptation scheme in the β-δ domain(exemplary embodiment). This non-linear adaptation scheme in the β-δdomain is shown in FIG. 2 along with a linear adaptation scheme. Fromequation (6), we can see that β attains its maximum value atδ_(m)=(Δ/2).

The maximum value is β_(m)=(Δ/4). Clearly, to avoid blocking artifacts,we must choose small Δ. This limits the capability of the noisesuppression.

In accordance with another aspect of the present invention, a linearadaptation in the β-δ domain can be written as:

$\begin{matrix}{\beta = \{ \begin{matrix}{{{\kappa\delta}^{\prime} + \chi},} & {{\delta^{\prime} = {\max\{ {\delta,\delta_{0}} \}}},} & {\delta < \delta_{1}} \\{0,} & \; & {\delta \geq \delta_{1}}\end{matrix} } & (7)\end{matrix}$where κ and χ are constants less than one. Thus, equation (7) representsa linear calculation of the discrimination value for the filter, whichdiscrimination value is then added or subtracted to/from the input frameto obtain the output frame.

In accordance with another aspect of the present invention, β should betruncated down to an integer to avoid trailing or blocking artifacts. Ingeneral, because the constants κ and χ are fixed numbers, equation (7)has to be implemented with 16-bit precision.

From equation (4), we can see a linear adaptation in the β-δ domain istranslated into a non-linear adaptation in the α-δ domain. Thenon-linear adaptation in the α-δ domain is shown in FIG. 3 along withthe optimal adaptation. From FIG. 3, we can see the non-linearadaptation closely approximates the optimal adaptation. Also we can see,the implementation of the non-linear adaptation in the α-δ domain israther complex as compared to the implementation of the same adaptationin the β-δ domain, in which the adaptation becomes linear.

In accordance with still another aspect of the present invention, fromequation (7), filters with different strength can be implemented byselecting different slopes. This is another advantage over thenon-linear adaptation of equation (6). From equation (6), the slope, atδ=0, is actually fixed regardless of the strength of the filter. Astronger filter can be implemented either by increasing the slope or byincreasing the maximum discrimination β. Equation (7) allows us toincrease one without changing the other. This allows the adaptationscheme be more manageable.

In accordance with yet another aspect of the present invention, fromequation (7), κ and χ are used primarily to adjust the filter strengthwhere δ is small, while δ₀ and δ₁ are used to adjust the filter strengthwhere δ is relatively large. δ₀ is also used to adjust the maximumdiscrimination value to allow as much noise reduction as possiblewithout disrupting the video image. Equation (7) allows us to adjust oneaspect of the filter without adversely affecting other aspects of thefilter. This allows the adaptation scheme be more manageable.

An example of linear adaptation in the β-δ domain is given as:

$\begin{matrix}{\beta = \{ \begin{matrix}{{( {\delta^{\prime} + 3} )\text{/}4},} & {{\delta^{\prime} = {\max\{ {\delta,5} \}}},} & {\delta < 7} \\{0,} & \; & {\delta \geq 7}\end{matrix} } & (8)\end{matrix}$where we choose κ=¼, χ=¾, δ₀=5, and δ₁=7.

From equation (8), we can see, the implementation of the linearadaptation, which maps to a non-linear adaptation in the α-δ domain, canbe performed with 8-bit precision. In other words, the exemplaryembodiment of the non-linear adaptation in the α-δ domain not onlyclosely approximates the optimal adaptation, but also reduces thecomputational complexity (note that the division can be implemented byright shifting).

The optimal adaptation scheme is established experimentally to allowmaximum noise suppression without adversely affecting the video content.A new implementation of MCTF is developed, which employs non-linearadaptation scheme to closely approximate the optimal scheme. The newimplementation also reduces the computational complexity if parametersare selected appropriately.

Turning to FIG. 5, shown therein is an exemplary embodiment of a methodfor filtering a video signal according to still another aspect of thepresent invention.

In element 51, the next frame in the sequence of frames is received bythe processor. We shall call this frame the current frame.

In element 52, the magnitude difference between the current frame and aprior filtered frame along the motion trajectories in the sequence offrames is determined.

In element 53, a discrimination value is determined based on themagnitude difference. This discrimination value can be obtained fromeither equation (6) or equation (7), although equation (7) isconsiderably easier to implement and manage and it also allows morenoise reduction without adversely affecting the video content.

In element 54, the current frame is modified using the discriminationvalue calculated above. The discrimination value is added to orsubtracted from the current frame to obtain the filtered version of thecurrent frame.

As the discrimination value is modified based on the magnitudedifference between the current frame and the immediately preceding framebefore each frame is “filtered,” the discrimination value thereby adaptsto the magnitude difference between frames.

In element 55, the strength of the filter is adjusted by modifyingcertain variables in the determination of the discrimination value.Various techniques can be employed for adjusting the strength of thefilter, such as: (1) by modifying the slope of the equation calculatingthe discrimination value β of the filter; (2) by modifying the maximumdiscrimination β; (3) in cases in which δ is small, by modifying theoffset of the equation χ; (4) in cases where δ is small by modifying theslope of the equation κ; (5) in cases where δ is relatively large bymodifying the cutoff value of the equation δ₁; and (6) in cases where δis relatively large by modifying the saturation value of the equationδ₀. These are only some examples of how to adjust the strength of thefilter.

In step 56, the maximum discrimination value is regulated by modifyingthe saturation value of the equation (e.g., δ₀ in some embodiments) toavoid trailing or blocking artifacts.

Although various embodiments are specifically illustrated and describedherein, it will be appreciated that modifications and variations of theinvention are covered by the above teachings and are within the purviewof the appended claims without departing from the spirit and intendedscope of the invention. For example, certain filter values are depicted,however, others are possible depending upon the application withoutdeparting from the scope of the present invention. Furthermore, theseexamples should not be interpreted to limit the modifications andvariations of the invention covered by the claims but are merelyillustrative of possible variations.

1. A method for filtering noise in a video encoding process comprising:determining a discrimination value of a filter based on a magnitudedifference between a current frame and a prior frame of a sequence offrames; and modifying the current frame with a noise reduction processorusing the determined discrimination value; wherein the current frame ismodified as follows:g(i,j,k ) =f(i,j,k )+β in which: f (i,j,k ) is the current frame inputto the filter; g (i,j,k ) is the current frame output from the filter;and β is the discrimination value of the filter; and wherein thediscrimination value is determined as follows:$\beta = \{ \begin{matrix}{{\delta - {( {1\text{/}\Delta} )\delta^{2}}},} & {\delta < \Delta} \\{0,} & {\delta \geq \Delta}\end{matrix} $ in which δ is the magnitude difference between thecurrent frame and the prior frame along motion traiectories; and Δ is apredetermined filter strength.
 2. The method according to claim 1,wherein the discrimination value is truncated to an integer prior tomodifying the current frame.
 3. The method according to claim 1, furthercomprising: adjusting the strength of the filter by modifying the slopeof the equation calculating the discrimination value β of the filter. 4.The method according to claim 3, further comprising: adjusting thestrength of the filter by modifying the maximum discrimination β.
 5. Anapparatus for encoding a video signal comprising: a first encoderreceiving the video signal; a second encoder receiving a plurality ofparameters from the first encoder to use in encoding the video signal;and a noise reduction preprocessor: receiving pixel data from the firstencoder; filtering the pixel data; determining a discrimination value ofa filter based on a magnitude difference between a current frame and aprior frame of a sequence of frames; modifying the current frame usingthe modified discrimination value; and providing preprocessed video tothe second encoder; wherein the noise reduction preprocessor modifiesthe current frame as follows:g (i,j,k ) =f (i,j,k )=β in which: f (i,j,k ) is the current frame inputto the filter; g (i,j,k ) is the current frame output from the filter;and β is the discrimination value of the filter; and wherein the noisereduction preprocessor determines the discrimination value as follows:$\beta = \{ \begin{matrix}{{\delta - {( {1\text{/}\Delta} )\delta^{2}}},} & {\delta < \Delta} \\{0,} & {\delta \geq \Delta}\end{matrix} $ in which δ is the magnitude difference between thecurrent frame and the prior frame along motion traiectories; and Δ is apredetermined filter strength.